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Abstract 

A fractional diffusion equation based on Riemann-Liouville fractional deriva- 
tives is solved exactly. The initial values are given as fractional integrals. 
The solution is obtained in terms of //-functions. It differs from the known 
solution of fractional diffusion equations based on fractional integrals. The 
solution of fractional diffusion based on a Riemann-Liouville fractional time 
derivative does not admit a probabilistic interpretation in contrast with frac- 
tional diffusion based on fractional integrals. While the fractional initial value 
problem is well defined and the solution finite at all times its values for t — > 
are divergent. 
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Anomalous subdiffusive transport appears to be a universal experimental phenomenon 
Examples occur in widely different systems ranging from amorphous semiconductors 
through polymers |§-|8| and composite heterogeneous films || to porous media |]I0,IT 



Theoretical investigations into anomalous diffusion and continuous time random walks have 
been a major focus of H. Schers research for many years [|T^,|T^,|,|r^] . The purpose of this 
paper is to discuss a theoretical approach based on the replacement of the time derivative 
in the diffusion equation with a derivative of noninteger order (fractional derivative). 

Many investigators have proposed the use of fractional time derivatives for subdiffusive 
transport on a purely mathematical or heuristic basis |l^-[T9|fl|20| . From the perspective 



of theoretical physics this proposal touches upon fundamental principles such as locality, 
irreversibility and invariance under time translations because fractional derivatives are non- 



local operators that are not invariant under time reversal [pl]| . These issues are generally 
avoided in heuristic and mathematical proposals, but were discussed recently in the context 
of long time limits and coarse graining It was found that fractional time derivatives 
with orders between and 1 may generally appear as infinitesimal generators of a coarse 
grained macroscopic time evolution PD[-^| . 

Differential equations involving fractional derivatives raise a second basic problem, re- 
lated to the first, that will be the focus of this paper. The second problem is whether to 
replace the integer order derivative by a Riemann-Liouville, by a Weyl, by a Riesz, by a 
Grunwald or by a Marchaud fractional derivative (see [25[-p8|j for definitions of these differ- 



ent derivatives). Different authors have introduced different derivatives depending on the 
physical situation pT| , p9| -pl|] . 

Given the basic objective of introducing fractional derivatives into the diffusion equation 
the present paper will be concerned with the equation 

D« + /(r, t) = C a Af(r, t) (1) 

where f(r, t) denotes the unknown field and C a denotes the fractional diffusion constant with 
dimensions [cm/s a ]. The fractional derivative operator, denoted as Dq + , is the Riemann- 
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Liouville derivative of order a and with lower limit a£l It is defined as [|25|,|28 

d 



dx 



(2) 



where 



r (a) 



(x-yy-'fiy) dy 



(3) 



is the Riemann- Liouville fractional integral with order a and lower limit a. Although many 
authors have investigated fractional diffusion problems it seems that equa- 

tion ([I]) has not been solved previously. In fact it was recently questioned whether an ap- 
proach using eq. ([]]) is consistent pi,29fl. It is the purpose of this paper to solve eq. (|Tj) 
exactly thereby establishing its consistency for appropriate initial conditions. 

Let me emphasize that eq. (|l|) differs from the popular equation introduced and solved 



m 



16fl . The latter equation is obtained by first rewriting the diffusion equation in integral 



form as 



f(r,t) = f 5(r) + C l [ A/(r,0 dt' 
Jo 



(4) 



where C% is the usual diffusion constant, 8(r) is the Dirac measure at the origin, and where 
the initial condition f(r, 0) = fo5{r) has been incorporated. Then the integral on the right 
hand side is replaced by a fractional Riemann-Liouville integral to arrive at the fractional 
integral form 

f(r,t) = f 6(r) + fit - O^A/Cr.O dt' = f 5(r) + C a (I% + Af)(r,t) (5a) 

r(") Jo 



or, upon differentiating both sides, at 





dt 



(5b) 



where C a is again a fractional diffusion constant. For a = 1 this reduces to eq. (||). The 
exact solution of eq. ([Sp is known and given by eq. ( |22|) below. 

In P5] , p6| it was shown that eqs. (^) have a rigorous relation with continuous time random 
walks of the kind investigated frequently by Harvey Scher |T^,0,||,|lOl . More precisely, eq. 



(H) was found to correspond exactly to a continuous time random walk with the long tailed 
waiting time density 



1 / t\ a -' ( t a 



where r is a time constant. Here E a ^(x) denotes the generalized Mittag-Leffler function 
defined by 

CO 

for all a > and b G C. For a = 1 this reduces to an exponential waiting time density. 
For < a < 1 these waiting time densities have a long tail decaying as ip(t) ~ t~ 1-a 
for £ — ► oo. Interestingly ip(t) ~ t a_1 diverges algebraically for t — > 0. It follows from 



refs. P^ , |^ , |^ . [24| . |21| that among the the waiting time densities with long tails the densities 
ip(t; a, To) represent important universality classes for continuous time random walks. 

Note that eqs. (|J) and ([I]) are not equivalent. The difference between eqs. (||) and (HD 
has to do with the initial conditions. An appropriate inital condition is found by analysing 
the stationary case. One finds that the fractional integral 

lS- a) f(r,0+) = f , a 6(r) (8) 

is preserved during the time evolution. This is a nonlocal initial condition. It implies the 
divergence of f(r,t) as t — > 0, as is characterstic for fractional stationarity p0|j2^] . 

Equation (|I|) with initial condition (|§p can be solved exactly by Fourier-Laplace tech- 
niques. Let the Fourier transformation be defined as 

*{f(r)}(q)= [ e l " r f(r)dr. (9) 

Fourier and Laplace transformation of eq. ([!]) now yield 

'<«■"> = dr*- (10) 

Inverting the Laplace transform gives 

f{q,t) = f , a t^E a ^-C a q 2 t a ). (11) 
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Setting q = shows that f(r,t) cannot be a probability density because its normalization 
would depend on t. Hence eq. (^]) does not admit a probabilistic interpretation contrary to 
eq. (§). 

To obtain f(r,t) it is advantageous to first invert the Fourier transform in eq. (|T0| ) and 
only later the Laplace transform. The Fourier transform may be inverted by noting the 



formula 37 



(27T) 



-d/2 I e iq-r 



m 



l-(d/2) 



K( d _ 2)/2 {m\r\) dr 



q 2 + m 2 



(12) 



which leads to 

f(r,u) = f 0!a (2nC a )- d / 2 



l-(d/2) 



ru 



a/2 



^a / V V 

with r = \r\. To invert the Laplace transform it is convenient to use the relation 

M {£{/(*)} («)}(!-*) 



(13) 



■M {/(*)}(*) 



r(i- s ) 



(14) 



between the Laplace transform and the Mellin transform 



M{f(t)}(s)= I t s -'f(t)dt. 

10 



(15) 



of a function /(£). Setting A = r/y/C^, X = a/2, v = [d — 2)/2 and /i = a(d — 2)/4 and 
using the general relation 



leads to 



M {x q g(bx p )} (s) = -b-^/ p g ( 

V V V 



\v > o) 



(16) 



M {f(r, u)} (s) = J -^{2nC a )- d ' 2 A 1 -^A~^l x M {K u {u)} {{s + /i)/A) . (17) 



A 



The Mellin transform of the Bessel function reads |38 



m {k v {x)} (s) = r- 2 r ( i±^] r ( s v 



(18) 
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Inserting this, using eq.(|14D, and restoring the original variables then yields 



M{f(r,t)}(s) 



a(r 2 7r) d / 2 \2^/C, 



T(l-s) 



(19) 



for the Mellin transform of /. Comparing this with the definition of the general if-function 



in eqs. (|23|) , (T2~4|) allows to identify the i7-function parameters as m = 0,n = 2,p = 2,q = 1, 
A 1 = A 2 = l/a, at = 1 - (d/2) + (1 - (1/a)), a 2 = (1 - (1/ar)), h = and Bx = 1 if 
(ad/2) — (a — 1) > 0. Then the result becomes 

2(l-(l/a)) 



0,Q 



tt02 
H 21 



a{r 2 7c) d / 2 \2^ c 
1 - 



(0,1) 



(20) 



This may be simplified using eqs. (pop, fl36|) and (^7]) to become finally 



/(r,t) = ; ?'° \„„ g- 20 



(r 2 7r) rf / 2 
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(a, a) 



\ 



(21) 



y - (d/2, l), (i,i) J 

This result should be compared with the known solution 

/ 

f(r,t) 



fo 



{r 2 7i) d / 



20 
2^12 



-H 



AC a t c 



\ 



(l,a) 
[d/2, 1) 



(22) 



1,1) 



of eq. (H) in which case f{r,t) is also a probability density. 

In summary this paper has shown that fractional diffusion based on Riemann-Liouville 
derivatives requires a fractional inital condition given by eq. (|8|). With this initial condition 
the fractional Cauchy problem can be solved exactly in terms of if-functions, and the solution 
is similar to the exact solution of the fractional integral form in eq. (Q. However contrary 
to eq. (|5|) whose solution is a probability density, and which is related to continuous time 
random walks, the solution of eq. (|I|) does not admit a probabilistic interpretation. 

The reader may ask why it is important for theoretical physics and chemistry to inves- 
tigate different forms of fractional diffusion. An answer was given already in p0|-f24|1, and 



recently again in |39| , |40[| . In these works it was found that fractional time derivatives arise 
generally as infinitesimal generators of the time evolution when taking a long-time scaling 
limit. Hence the importance of investigating fractional equations arises from the necessity 
to sharpen the concepts of equilibrium, stationary states and time evolution in the long time 
limit. 



APPENDIX: H-FUNCTIONS 



The if -function of order (m,n,p,q) G N 4 and parameters Ai G R+(z = 1, . . . ,p), Bi G 
(i — 1, . . . , q), di G C(i = 1, . . . ,p), and bi G C(i = 1, . . . , q) is defined for z G C, z ^ 



by the contour integral [fjd| -|45 

/ 

z 



Trm,n 



\ 



(ax,Ai), ... , (a p ,A p ) 



(fti.SO,... ,(b q ,B qJ 



2ni 



T](s)z s ds 



(23) 



where the integrand is 



Y[r(bi + B iS ) jjr(i-ai-As) 



?7(s) 



i=l 



i=l 



r(a i + A jS ) 11 r(i -bi-Bis) 



(24) 



i=n+l 



i=m+l 



In (23) z s = exp{— slog \z\ — i&rgz} and argz is not necessarily the principal value. The 
integers m,n,p,q must satisfy 



0<m<g, < n < p 



(25) 



and empty products are interpreted as being unity. The parameters are restricted by the 
condition 



where 



P a n F b = 



P a = {poles of T(l - en - A iS )} = <J - — eC:i' = l,...,n;jfcgN 



(26) 



As 



¥ b = {poles of r(6 < + B iS )} 



-h - k 
Bi 



G C : i = 1, . . . ,m;fcGN 



(27) 



are the poles of the numerator in (24). The integral converges if one of the following condi- 
tions holds 



£ = £(c-ioo,c + ioo;P a ,P 6 ); | arg^| < C+/2; C > (28a) 
£ = £(c-ioo,c + zoo;P a ,P fe ); \ eagz\ = C+/2; C > 0; cD < — Re F (28b) 



£ = £(-oo + 171, -00 + 172; P a , P&); -D > 0; 0<|^|<oo (29a) 
£ = £(-oo + i7i,-oo+i7 2 ;P a) P 6 ); D = 0; < \z\ < E' 1 (29b) 
£ = £(-oo + i7i,-oo + i7 2 ;P tt ,P 6 ); .D = 0; \z\ = E' 1 ; C > 0; ReF < (29c) 



£ = £(00 + 271,00 + 172; P a , D"b); £> < 0; < |z| < 00 (30a) 
£ = £(oo+i7i,oo + i72;P ,P 6 ); D = 0; |*| > F" 1 (30b) 
£ = £(oo+i7 1 ,oo + i7 2 ;P a ,P 6 ); D = 0; |z| = F -1 ; C > 0; ReF < (30c) 

where 71 < 72. Here C(zi, z 2 \ Gi, G 2 ) denotes a contour in the complex plane starting at Z\ 
and ending at z 2 and separating the points in Gi from those in G2, and the notation 

n p m q 



1=1 
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i=n+l 
V 



i=l 



i=m+l 



i=l i=l 

1=1 8=1 
q P 



(32) 

(33) 

(34) 
i=i i=i 

was employed. The F-functions are analytic for z 7^ and multivalued (single valued on 



the Riemann surface of log z). 

A change of variables in (|23"D shows 
/ 



Trm,n 



V 



, (cip, Ap) 



(ai, Ai), 

(b 1 ,B l ),... ,(b q ,B q ) J 



\ 


( 

1 


(1 


-61, -B0,. 


•,(1 


TTn,m 

~ n q,p 

) 










V 


(1 


- oi, AO, • 


..,(1 



(35) 



which allows to transform an ^-function with D > and aigz to one with D < and 
arg(l/z). For 7 > 
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Tjm,n 
H p,Q 



1 (ai,Ai),... ,(a p ,A p ) ^ 

z 

\ (b 1 ,B 1 ),... ,(b q ,B q ) j 



( 



Tjm,n 



(01,7^1),... ,(a p ,7A p ) ^ 
(61,750,... ,(b q ,-fB q ) J 



(36) 



while for 7 e 



7 rrm,n 



V 



(ai,Ai),... , (a p ,A p ) ^ 
(b 1 ,B 1 ),...,(b q ,B q ) j 



Tjm,n 



(ai + 7^1, Ax), . . . , (ap + 7^ P , 40 ^ 
(h + iB^Bi),... ,(b q + >yB q ,B q ) j 



(37) 



holds. 
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